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ABSTRACT: With the large-scale integration of renewable 
energy and power electronic devices, system voltage support 
strength decreases, thus threatening the system's safe and stable 
operation. In the homogeneous scenario where power 
electronic devices integrated into the system have the same 
dynamics, a theoretically rigorous and highly stable 
quantitative method for quantifying system voltage support 
strength can be formed based on the generalized short-circuit 
ratio(gSCR) and the device critical short-circuit ratio(SCRo) In 
the heterogeneous scenario where power electronic devices 
integrated into the system have weakly different dynamics, 
system voltage support strength can be quantified by the 
first-order approximation of gSCR and SCRo, based on some 
special dynamic characteristics of devices and power grid. 
However, there is a lack of unified calculation principles at the 
theoretical level. To this end, this paper focuses on quantifying 
system voltage support strength under small-signal stability and 
discovers that the multi-infeed system can be approximately 
decoupled into multiple low-dimensional systems. On this 
basis, this paper proposes the concept and approximate 
calculation method of eigen-subsystems and interprets their 
Based on the 


eigen-subsystems, the general calculation principle of gSCR 


physical significance. concept of 
and SCRo were unified, which achieves strength quantification 
of heterogeneous multi-infeed systems. Additionally, the 
specific calculation methods for gSCR and SCRo are provided 
in typical scenarios, such as grid-following converters under 
non-rated operating conditions, reverse active power output, 
and considering grid-forming devices. Finally, the effectiveness 
of the calculation principles and methods is verified in several 
Cases. 
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摘要 : 随 着 新 能 源 及 电力 电子 设备 的 大 规模 接 入 ,系统 电 
支撑 强度 降低 , 新 型 电力 系统 的 安全 稳定 风险 增加 。 在 局 
电力 电子 设备 接 入 场景 下 , 利用 电网 广义 短路 比 和 设备 临 
短路 比 可 以 形成 理论 严谨 、 与 稳定 强 相 关 的 系统 电压 支撑 
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度量 化 方法 ; 在 弱 异 构 场 景 下 , 基于 设备 和 电网 一 些 特殊 动 
态 特性 , 利用 广义 短路 比 和 设备 临界 短路 比 的 一 阶 近似 可 以 
实现 强度 量化 , 但 缺乏 通用 的 计算 原理 。 为 此 ， 本 文 聚焦 小 
扰动 稳定 维度 下 电压 支撑 强度 的 量化 问题 , 首先 利用 电力 电 
子 多 馈 入 系统 可 近似 解 耦 为 多 个 低 维系 统 的 规律 , 提出 特征 
子 系统 的 概念 和 近似 计算 方法 , 并 诠释 J ? 
基于 特征 子 系统 , 提出 了 广义 短路 比 及 其 临 
计算 原理 ,实现 了 异 构 多 馈 入 系统 的 强度 量化 。 
跟 网 型 电力 电子 设备 在 非 额 定 运行 点 、 有 功 功率 反 向 以 及 包 
含 构 网 型 电力 电子 设备 等 几 种 典型 场景 ,给 出 了 具体 计算 公 
式 。 最 后 ， 算 例 验证 了 所 提 原 理 和 方法 的 正确 性 和 有 效 性 。 


关键 词 : 广义 短路 比 ， 设 备 临 界 短路 比 ， 特 征 子 系 统 
0 引言 


随 着 新 能 源 的 快速 发 展 ， 电 力 系统 逐步 演化 成 
含有 高 比例 新 能 源 及 电力 电子 设备 的 新 型 电力 系 
统 中 。 风 电 、 光 伏 等 新 能 源 集中 并 网 容量 大 、 与 
主 网 电 / 距离 远 ， 上 且 常 采用 跟 网 型 控制 策略 而 对 电 
网 支撑 强度 弱 ， 使 和 租 系 统 电 压 支 撑 强 度 变 低 ， 受 扰 
后 安全 稳定 风险 增加 BH， 吸 需 新 型 电力 系统 电压 
支撑 强度 的 科学 量化 和 计算 方法 。 其 中 ， 考 虑 到 低 
惯量 特征 的 系统 强度 问题 可 以 分 开 考虑 ， 后 文 为 表 
述 方便 ， 将 电压 支撑 强度 简写 为 “系统 强度 ”。 

大 规模 电力 电子 设备 接 入 后 ， 系 统 强度 的 分 析 
计算 存在 以 下 难点 : 实际 系统 中 设备 及 中 间 节 点 等 
数量 较 多 ， 且 设备 与 网 络 动态 交互 复杂 ， 高 维 动态 
特征 难以 简化 分 析 ; 新 型 电力 系统 中 设备 种 类 众多 ， 
如 控制 策略 和 参数 等 各 异 的 跟 网 型 变 流 器 、 构 网 型 
变 流 器 或 调 相 机 等 I0 习 ， 显著 增 强 系统 的 复杂 程度 ; 


T 
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新 能 源 运行 方式 多 样 ， 需 要 考虑 实际 工 况 ， 也 给 系 
统 强度 量化 带 来 巨大 挑战 。 
为 分 析 系 统 强 度 ， 国 内 外 多 个 权威 机 构 和 组 织 


都 默认 电网 强度 与 计 及 设备 动态 的 系统 强度 之 间 
存在 正 相 关 性 ， 从 而 可 找到 描述 系统 强度 的 指标 ， 

例如 : CIGRE 多 馈 入 短路 比 (multi-infeed SCR， 
MSCR )037 0605 、 位 置 相关 短路 比 (site-dependent short 


SDSCRJ)JH 、 多 场 站 短路 比 
(multi-renewable energy station Short circuit ratio ， 

MRSCR)N 7、 等 效 短路 比 (equivalent circuit-based 
SCR ，ESCR)09 和 带 交 互 因 子 的 短路 比 (the SCR 


circuit ratio, 


with interaction factors, SCRIEF)GI 等 指标 。 然 而 ， 上 
述 指标 均 基 于 系统 静态 方程 ， 与 稳定 性 关联 弱 ， 


物理 意义 不 清晰 , 即使 量化 设备 外 特性 相似 的 电力 
电子 多 馈 入 系统 强度 也 是 如 此 ; 此 外 ， 上 述 指标 基 
于 短路 电流 ， 而 有 些 电 力 电 子 设备 可 以 不 提供 短路 
电流 ， 但 它们 对 稳定 性 影响 不 容 忽 略 ， 因 此 ， 传 统 
系统 强度 量化 方法 在 复杂 场景 下 的 误差 很 大 。 
进一步 ， 文 献 [7]、[18]-[20] 综 合 考 虑 电网 强度 
与 设备 动态 以 量化 系统 强度 ， 基 于 系统 的 灵敏 度 方 
程 提 出 了 广义 短路 比 (generalized short circuit ratio， 
gSCR) 指 标 。 广 义 短路 比 指标 与 系统 强度 存在 解析 
关系 ， 与 稳定 性 强 相 关 。 在 同 构 系 统 〈 接 入 设备 的 
动态 特性 相似 ) 中 ， 多 电力 电子 设备 馈 入 系统 〈 简 
称 “ 多 馈 入 系统 ”) 的 整体 动态 可 解 耦 为 多 个 等 效 
单 馈 入 系统 的 动态 ， 基 于 等 效 单 馈 入 系统 ， 可 得 到 
电网 广义 短路 比 和 设备 临界 短路 比 ， 实 现 原 系统 的 
强度 量化 ; 在 弱 异 构 场景 下 ， 如 外 特 a 
直 驱 风机 器、 光伏 加、 构 网 型 设备 L000N04、 储 能 系 
统 接 入 PH 或 设备 在 非 额定 运行 点 [ 人 可 以 将 
多 样 化 设备 动态 看 作对 同 构 系 统 中 设备 动态 的 摄 
动 或 网 络 动态 的 修正 ， 从 而 修正 广义 短路 比 和 设备 
临界 短路 比 实现 强度 量化 。 上 述 计算 关键 为 求解 能 
够 近似 原 系统 稳定 性 的 低 维 单 馈 入 系统 ， 在 此 基础 
上 实现 gSCR 指标 及 临界 值 的 定义 和 计算 。 然 而 ， 
现 有 研究 给 出 的 计算 方法 仅 适用 于 固定 的 场景 ， 划 
理论 缺乏 普 适 性 。 

为 此 ， 本 文 引 入 特征 子 系统 的 概念 ， 保 留 设 备 
与 网 络 的 动态 特性 , 统一 新 型 电力 系统 广义 短路 比 
及 系统 强度 的 计算 原理 。 首先 , 引入 了 特征 子 系统 
的 概念 并 提出 其 求解 方法 ; 其 次 ， 基 于 特征 子 系 
统 ， 提 出 了 广义 短路 比 与 设备 临界 短路 比 的 计生 
方法 ， 统 一 了 多 种 场景 下 的 基于 广义 短路 比 的 强 
度量 化 方法 ; 最后， 通过 算 例 验证 了 方法 的 有 效 
性 。 


1 预备 知识 
1.1 矩阵 的 特征 子 空 

2nx2n 矩阵 Who 的 特征 值 、 
特征 和 矩阵 定义 如 下 ; 

定义 1.1031 (特征 值 及 特征 向 量 ): 矩阵 的 特征 
值 为 特征 方程 dt(M 一 4L,) 必 的 根 ( 式 中 4 为 27 
维 单位 矩阵 ); x 和 ?分别 为 矩阵 与 特征 值 / 对 应 的 
右 和 左 特征 向 量 ， 满 足 : 


Wl 


s 间 与 广义 特征 矩 阵 
特征 子 空 间 及 广义 


Mx=Ax,y” M,=Ay" 人 
其 中 ,x 和 了 均 为 27 维 列 向 量 。 
当 Mo 的 元 素 为 s 的 函数 时 ，?ho 的 特征 值 为 
的 函数 ， 拓 展 为 “特征 函数 ” 
定义 1.254《〈 特 征 子 空间 和 广义 特征 矩阵 ): 令 
R(X,)= {Xou:ue Cc”™"} 、R(Y,)= {You:u e c”™")} 分 别 
为 矩阵 Mp 的 右 、 左 特征 子 空间 ， 分 别 由 Xo、 的 
列 向 量 张 成 ， 其 中 ， 不 失 一 般 性 ， 归 一 化 后 满足 
YoYXo=1; 与 特征 子 空间 对 应 ， 存 在 唯一 的 矩阵 Lo， 
Lo 的 特征 值 是 No 的 特征 值 ， 称 之 为 “广义 特征 拖 
阵 ”， 满足 下 式 。 


MX = XLo, Ye Mo, = Ly (2) 
进一步 ，Zo 可 由 下 式 求 解 得 到 ; 
L =Y" MoX, (3) 


值得 一 提 的 是 , 广义 特征 和 矩阵 及 特征 子 空间 分 
别 是 特征 值 及 特征 向 量 的 推广 ， 满 足 : 若 ("， 力 是 
广义 特征 矩阵 Zo 的 特征 对 ( 式 中 ”为 Zo 的 右 特征 向 
量 )， 则 (CXov， 力 为 矩阵 Mo 的 特征 对 ， 其 中 ， 特 征 
值 是 维 数 为 1 的 特殊 广义 特征 矩阵 ， 特 征 向 量 张 
成 的 空间 是 特征 子 空间 的 特例 , 对 应 Lo 为 对 角 阵 。 
进一步 ， 由 上 述 分 析 可 知 ， 低 维 广义 特征 和 矩阵 
Lo 的 特征 值 为 矩阵 Mo 的 特征 值 ， 则 Lo 可 理解 为 通 
过 数学 变换 得 到 的 Mo 的 解 耦 子 系统 。 为 方便 理解 
及 后 文 推 导 ， 下 文 举例 说 明 。 
假设 矩阵 Mo。 有 一 组 特征 子 空间 RX )、(Y) 
(i=1,...,n), 其 中 ，X 二 0 Xi)、 Y=(ya, p12), 对 应 n 
个 二 维 广义 特征 矩阵 L;， 如 下 所 示 。 则 令 咎 CO .…. 
XK) 和 Y= ... 功 )， 可 将 11o 转换 为 矩阵 卫 ， 其 对 角 
块 为 n 个 二 维 广 义 特 征 和 矩阵 L;:， 表 示 为 : 
L=Y"MX=diag(L) (4) 
式 中 : 二 = 了 ”MMX 
由 于 X、 工 的 列 向 量 都 线性 无 关 ， 和 矩阵 工 与 
Mo 的 特征 值 完 全 相同 , 其 中 工 ;的 特征 值 为 Wo 的 部 
分 特征 值 ， 则 L; 可 理解 为 数学 变换 得 到 的 Mo 的 n 
个 解 起 子 系统 ， 分 别 利用 R( 关 )、R(Y) 即 可 求 得 。 
特别 的 ,车 Mo 及 n 个 二 重 特征 值 4( 二 1,…,n) 且 存在 
两 个 线性 无 关 的 特征 ， 其 张 成 的 空间 恰好 为 特 
征 子 空间 ， 此 时 计算 可 得 广义 特征 矩阵 L 二 4。 
1.2 ”特征 子 空间 摄 
考虑 如 下 受 扰 和 矩阵 : 
M=M,+AM (5) 
其 中 ， 和 矩阵 AM 为 矩阵 Mp 所 受 扰 动 。 
若 扰动 矩阵 AM 足够 小 ， 克 的 特征 子 空 间 与 
Mo 的 特征 子 空间 近似 相同 叫 。 根 据 该 性 质 ， 可 求 
得 如 下 扰动 后 让 广义 特征 德 阵 的 一 阶 近似 P4。 


汪 
这 
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句 量 - 


第 期 刘 晨 曦 等 ， 基于 特征 子 


引 理 1 p4 (特征 子 空间 摄 动 ): 矩阵 特征 子 空 
闻 的 零 阶 近似 为 只 羽 ) 、 噶 ( 肪 ) ， 六 广义 特征 矩阵 
态 的 一 阶 近似 可 由 下 式 求解 得 到 。 


L, = Y” MX, = L,+(Y AMX,) 


其 中 ，Zo 为 而 的 零 阶 近似 。 

由 引 理 1 可 知 ， ?0 的 特征 子 空间 可 以 逼近 万 
的 特征 子 空 间 。 扰 动 越 小 ， 两 个 矩阵 的 特征 子 空间 
越 接近 ， 二 广义 特征 矩阵 一 阶 近 似 的 误差 越 小 。 最 
然 ， 当 友 的 特征 子 空间 计算 非常 复杂 时 ， 可 通过 构 
造 合 适 的 Mo 简化 计算 。 为 方便 理解 ， 以 两 维 矩 阵 
为 例 说 明 上 述 结论 ， 见 附录 A。 
1.3 We。 


(6) 


定 动态 系统 4, 其 闭环 传递 浮 数 矩 阵 为 六 (9)， 


令 严 =T(s) (简称 “系统 矩阵 ”)， 则 
det( 故 )=0 为 特征 方程 ， 培 的 特征 函数 行列 式 能 够 
决定 4 的 稳定 性 ， 而 由 式 可 知 ， 克 可 解 看 为 多 个 
低 维 广义 特征 矩阵 , 每 个 广义 特征 矩阵 的 特征 函数 
都 是 方 的 特征 函数 。 因 此 ， 针 对 关注 的 动态 系统 部 
分 模 态 ， 存 在 对 应 特征 函数 与 广义 特征 矩阵 到 (k 
E {1,...,n)), 决定 了 动态 系统 在 关注 模 态 下 的 稳定 
性 ;与 五 对 应 ， 存 在 低 维 子 系统 44， 其 闭环 传递 


函数 为 系 。 记 A4i 为 动态 系统 4 的 特征 子 系统 ， 可 
近似 4 在 关注 模 态 下 的 稳定 性 ， 具 体 为 : 

定义 3.1 特征 子 系统 ): 设 动态 系统 4， 闭 环 
传递 函数 为 T!i(s)， 令 丰 =T(s)， 记 关注 模 态 下 应 的 
特征 子 空间 和 广义 特征 矩阵 分 别 为 (,)、 ( 评 ) 


和 工 , ， 满 足 式 ， 闭 环 传递 函数 为 4 的 低 维 子 系 
统 4k， 称 之 为 动态 系统 4 的 特征 子 系统 。 
L =Y MYX,,ke{l..,n) 
值得 注意 的 是 ， 特 征 子 系 统 4 能 够 近似 动态 
系统 4 在 关注 模 态 下 的 稳定 性 ; 此 外 ， 根 据 式 (7) 


可 知 荆 为 详解 艳 得 到 的 低 维 广义 特征 矩阵， 且 后 
文 将 说 明 特 征 子 系统 具有 实际 物理 意义 及 动态 特 
性 ， 因 此 特征 子 系统 4 的 物理 意义 为 动态 系统 4 
的 解 耦 子 系统 。 此 外 ， 由 于 系统 动态 特性 多 样 ， 系 
统 模 态 很 多 ， 对 应 动态 系统 有 多 个 特征 子 系统 ， 可 
分 别 近 似 不 同 模 态 下 的 稳定 性 。 
1.4 用 到 的 一 些 术语 
为 了 后 文 描述 ， 本 节 
语 说 明 。 
1) 


(7) 


结合 动态 系统 给 出 几 个 术 


系统 矩阵 (闭环 系统 ): 动态 系统 的 闭环 传 
递 函数 逆 矩阵 ， 其 行列 式 能 够 决定 动态 系 


2j 


系统 的 广义 短路 比 计 算 原 理 


统 稳定 性 ; 


广义 特征 矩阵 Lo: 


特征 函数 为 系统 矩阵 


(关注 模 态 下 ) 特 征 函 数 的 矩阵 称 为 广义 特 


征 和 矩阵 ， 


3) 特征 子 空 


下 ) 广 义 特 征 


空间 ， 满 
4) 
模 态 下 广 


维 数 较 低 ， 


足 式 -; 


中 广义 特征 矩阵 决 
定 了 动态 系统 在 关注 模 态 下 的 稳定 性 ; 


间 9( 5)、 


R( 叶 ) : 与 (关注 模 态 


原 系统 的 解 厢 子 系统 。 


5) 


同 构 ( 异 构 ) 系 统 : 
似 (存在 差异 ) 的 多 馈 入 系 
存在 差异 的 设备 包含 


和 矩阵 对 应 , 可 将 系统 矩阵 降 阶 
为 广义 特征 矩阵 的 列 向 量 集合 为 特征 


子 


特征 子 系统 : 闭环 传递 函数 逆 算 阵 为 关注 
义 特征 矩阵 的 低 维 
近似 动态 系统 在 关注 模 态 下 的 稳定 性 ; 具 
实际 物理 意义 与 动态 特性 ， 


子 系统 ， 能 够 


物理 意义 为 


电力 电子 设备 外 特性 相 
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1 所 示 电 力 电子 设备 接 入 的 新 型 电力 


系统 ， 包 含 ntm 个 设备 母线 ,r 个 中 间 无 源 母 线 和 


无 穷 大 母线 ( 若 有 构 网 


其 
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2.1 系统 动态 建 模 


首先 ， 考 虑 所 


在 全 局 xy 坐标 系 下 ， 
性 化 模型 ， 得 到 系统 的 闭环 特征 方程 。 
线路 阻抗 比 相同 ， 从 设备 端口 


往 交 流 电网 看 ， 在 
压 微 增 量 AUw(s) 和 
可 表示 为 3: 


全 局 同步 4 


设备 , 无 穷 大 母线 可 
1 设备 包含 直 驱 风机 、 光 伏 、 调 相机 和 储 能 系统 ; 
来 分 , 设备 可 以 分 为 基于 PLL 的 跟 网 型 
变 流 器 、 基 于 虚拟 同步 机 的 构 网 型 变 流 器 


建 并 图 


以 没有 )。 


[10][11][12] 


1 所 示 系 统 的 线 


标 系 下 的 电网 端口 电 
电流 微 增 量 ALw(s) 的 灵敏 度 方程 


AL, =(B, ®6(s)+ Be GAS))AU。 人 
Zu(5) 


B B 


Ll(ntm+r) 


工 11 
3 Oo 

B; = ,AS) = 2 | 

B … BB " 
L(n+m+r)l L(ntm+r)(ntm+r) 
es r+l s/(@, +z) 1 

re 

其 中 ，Bi 仪 包含 线路 感 纳 及 电导 信息 ; Bc 仅 包含 


线路 容纳 信息 ; t=R/X 为 线路 阻抗 比 A 表示 变量 
的 增 量 ， @ 表 示 Kronecker 积 ; wo 为 交流 电网 的 同 
步 转速 ，diag(:) 表 示 对 角 珑 阵 。 
其 次 ， 设 备 端 口 的 电压 微 增 量 ACo(9) 和 电流 微 
增 量 ALwy(s) 的 灵敏 度 方程 可 表示 为 L081: 

了 ma (s) 


人 
其 中 ，Yisri(s) ~ sn) 分别 为 n+7 台 设 备 的 导 
纳 传递 函数 矩阵 ，0>ER2” >， 元 素 均 为 0。 

最 后 ， 结 合式 和 ， 利 用 行列 式 可 表示 图 1 所 示 
系统 的 闭环 特征 方程 为 : 

det (T, (s)) = det(Yia(s)+Y,,(s))=0 (10) 

其 中 ，7Tiei(s) 为 闭环 传递 函数 逆 和 矩阵 。 
2.2 ”基于 广义 短路 比 的 系统 强度 量化 思路 回顾 
通过 求解 式 所 示 闭 环 特征 方程 的 根 可 量化 系 
统 稳定 裕 度 ， 然 而 ， 由 于 系统 中 设备 与 中 间 无 源 母 
线 数 目 较 多 导致 me(9) 维 数 较 高 ,计算 困难 。 为 此 ， 
考虑 寻找 近似 实际 系统 稳定 性 的 简单 系统 ， 在 此 基 
础 上 计算 广义 短路 比 ， 以 实现 系统 稳定 裕 度 分 析 ， 
从 而 实现 强度 量化 。 

针对 同 构 系 统 ， 多 馈 入 系统 的 整体 动态 可 解 耦 
为 多 个 等 效 单 馈 入 系统 的 动态 。 基 于 等 效 的 单 馈 入 
系统 ， 就 能 形成 融合 电网 广义 短路 比 gSCR 与 设备 
临界 短路 比 SCR。 的 系统 强度 和 稳定 性 判 据 ， 而 且 
具有 电网 动态 与 设备 动态 分 离 的 特点 。 具 体 为 : 

利用 式 - 所 示 gSCR 和 SCRo 度 量 系统 强度 ， 其 
中 : gSCR 仅 与 交流 电网 参数 和 设备 容量 参数 有 关 ， 
SCRo 与 设备 动态 特性 有 关 , 反映 设备 对 电网 短路 比 
的 耐 受 能 力 ， 它 也 是 系统 广义 短路 比 临界 值 
CgSCRo。 


1 1 
SSCR min 4 人 (SB ) =0” EC (1 1 ) 


其 中 ，SBg04m =diag(SB1,.., Sptm)，SB1.. SB(ntm) 分 别 
为 ntm 台 设 备 的 额定 容量 ; Be 为 忽略 电容 动态 且 
仅 保留 设备 母线 的 网 络 导 纳 矩阵 ; Zer Barl; 5() 
表示 求解 矩阵 的 最 大 奇异 值 。 
gSCR > CgSCR, = SCR, (12) 
根据 式 可 判断 系统 是 否 稳定 。 实 际 中 ， 为 应 对 
一 些 不 确定 性 因素 ， 可 要 求 广义 短路 比 不 仅 大 于 
SCRo， 还 需要 根据 现场 经 验 多 留 一 些 裕 度 。 因 此 ， 
根据 式 计算 gSCR 和 SCR。 指标 的 相对 值 p%6， 以 判 
断 系 统 是 否 满足 稳定 裕 度 要 求 。 例 如 ， 需 保留 至 少 
po% 的 裕 度 : 
gSCR— SCR, 


Bp% 
SCR, 

直 得 注意 的 是 ， 根 据 式 中 的 定义 可 知 ， 广 义 短 

路 比 与 阻抗 矩阵 的 最 大 奇异 值 有 关 ， 而 与 短路 容量 

没有 直接 关系 ， 它 本 质 是 反映 多 端口 交流 网 络 中 电 

压 扰动 对 电流 扰动 的 最 大 灵敏 度 D。 

如 果 非 同 构 系 统 ， 如 接 入 设备 包含 控制 策略 及 
参数 存在 差异 的 跟 网 型 变 流 器 、 构 网 型 变 流 器 、 调 
相机 等 ， 利 用 矩阵 摄 动 理论 对 gSCR 和 SCRo 进行 
近似 即 可 。 具 体 地 ， 根 据 摄 动 的 来 源 与 设备 电网 动 
态 分 离 特性 ， 可 分 为 两 种 : 1) 设备 动态 : 跟 网 型 
变 流 器 外 特性 存在 差异 时 ， 认 为 是 对 同 构 系 统 设备 
动态 的 摄 动 , 此 时 gSCR 的 计算 方法 不 变 , 而 SCR 
取决 于 设备 动态 特性 的 加 权 平 均 ， 详 细 可 参考 文献 
[7][19][21]; 2) 网 络 动态 : 当 考 虑 构 网 型 变 流 器 
或 调 相 机 接 入 等 场景 ， 可 认为 是 对 网 络 修正 ， 此 时 
SCRo 的 计算 方法 不 变 , 而 gSCR 需 考 虑 构 网 型 变 流 
器 和 调 相机 等 因素 ， 详 细 可 参考 文献 [9][11][12]。 
2.3 拟 解决 的 问题 

综 上 所 述 ，gSCR 指标 及 临界 值 的 定义 和 计算 ， 
其 关键 步骤 是 寻找 能 够 近似 系统 稳定 性 的 低 维 单 
馈 入 系统 ， 即 本 文 的 特征 子 系统 。 然 而 ， 现 有 研究 
仅 适 用 于 单一 典型 场景 下 的 问题 ， 方 法 缺乏 普 适 性 ， 
导致 对 广义 短路 比 的 理解 相对 复杂 。 为 此 ， 本 文 利 
用 更 具 物 理 意义 的 特征 子 系统 概念 ， 以 便于 理解 广 
义 短路 比 的 计算 原理 。 主 要 讨论 两 个 问题 : 

问题 1， 寻 找 近似 高 阶 系统 模 态 的 低 维 解 耦 子 
系统 ， 即 后 文 的 特征 子 系统 ， 其 理论 基础 是 什么 ? 

问题 2: 基于 特征 子 系统 ， 广 义 短路 比 及 其 
界 值 的 通用 计算 原理 是 什么 ? 


3 ”特征 子 系统 及 其 近似 方法 


结合 定义 3.1 可 知 , 对 给 定 动态 系统 及 其 闭环 
传递 函数 逆 矩 阵 ， 存 在 与 其 关注 模 态 对 应 的 低 维 
特征 子 系统 与 广义 特征 矩阵 。 甚 中， 特征 子 系统 


x100% > B,% (13) 


也 


Ew 


第 其 刘 晨 赚 等 ; 基于 特征 子 系统 的 广义 短路 比 计 算 原理 5 


的 闭环 传递 函数 为 对 应 广义 特征 矩阵 的 着 ， 可 近 Zou(s) 

似 动态 系统 在 关注 模 态 下 的 稳定 性 。 因此， 求解 式 中 ， Jo (可 = 本 

与 关注 模 态 对 应 的 实际 系统 闭环 传递 函数 逆 矩 阵 yo (s 

(简称 “实际 系统 矩阵 ”) 的 广义 特征 矩阵 (简称 de a 

“关注 模 态 下 广义 特征 矩阵 ”7， 可 得 到 特征 子 系 i 

统 。 进 一 步 ， 基 于 特征 子 系统 可 实现 实际 系统 的 BB mr (BB) 

强度 量化 ， 其 中 ， 特 征 子 系统 维 数 较 低 ， 分 析 和 a 

计算 都 更 加 简单 。 上 述 分 析 思路 的 关键 是 获得 特 Re 

征 子 系统 ， 即 求解 广义 特征 和 矩阵。 ER 
为 此 ,本 章 首先 基于 设备 特性 并 利用 和 矩阵 的 舍 L 。 


进一步 考虑 实际 工 况 ， 台 设 备 实 际 工 
况 与 其 动态 解 耘 。 有 具体 地 ， 用 包含 设备 实际 工 况 信 
息 的 矩阵 (C2/ 户 @ 天 左 乘 则 me 可 转化 
为 如 下 所 示 的 Tope(9)。 令 应 = T9146(s)，AMM 的 特 
征 函 数 决 定 了 图 1 所 示 系 统 的 稳定 性 。 


尔 补 变换 简化 实际 系统 和 矩阵 ， 其 次 ， 基 于 引 理 1 
获得 系统 矩阵 的 近似 特征 子 空间 ,并 得 到 对 应 的 低 
维特 征 子 系统 。 
3.1 “系统 矩阵 及 其 简化 

考虑 图 1 所 示 系 统 包含 了 ntm 台 设 备 ,其 中 m 


台 是 构 网 型 变 流 器 和 调 相机 等 ， 其 特点 是 设备 外 特 ot (VU/P) 1 )T (5) = eet(T() aa 
性 近似 为 含 内 电抗 的 电压 源 , 建 模 中 这 m 台 设 备 可 =det (Yi (s)+ Yr" (s))=0 


等 效 为 对 地 支 路 ， 等 效 导 纳 为 Br za (j=1,...,m); n ey fo 
台 设备 是 跟 网 型 新 能 源 变 流 器 与 常规 储 能 变 流 器 。 7 EA) 
等 设备 。 U0? U? . 
将 实际 系统 看 成 由 两 部 分 组 成 , -部 分 为 等 效 了 (*) ”B®6(s)+ pdiog (80)+A(Bc))® ps) 
网 络 ,包含 实际 网 络 与 m 台 构 网 设备 等 效 导 纳 ; 另 


一 部 分 包含 n 台 跟 网 设备 动态 模型 ， 表 示 为 ; i 0 
电压 、 输 深 。 
Ts (s) Yan(s) + els) 09 值得 注意 的 是 , 当 设备 均 在 额定 运行 点 , 设备 
| 端口 电压 均 为 1p.u.， 输 出 有 功 功率 分 别 等 于 自身 
T(s) T,(s) 额定 容量 Sgi(i=1,...,n)， 此 时 CPP diag(U2/P)) 
Yi (s)= (Bi +B, saa)®E(s)+ B® ns) =diag(SB71)=SBn1, 则 设备 在 额定 运行 点 是 实际 工 况 
0, 下 的 特例 。 
B, jpr = diag (Bi 8) 3.2 ”特征 子 系统 近似 方法 
0, 当 n 台 等 效 设备 动态 模型 同 构 时 ， 等 效 导 纳 甜 
Jan (s) 阵 (Z2/PB wu 8 卫 关 于 自身 二 重 特征 值 4(i=1,.….,n, 
yr (s)= 为 最 小 正 特征 信 ) 的 特征 子 空间 只 多 )、 中 区 为 
Ysn (5) 的 特征 子 空间 。 利 用 R(X,)、R(Y,) (ke {1,...n)) 
2(m+r) A | > 上 AR 
可 求解 实际 系统 的 多 个 特征 子 系统 ， 分 别 近似 实际 
其 中 ，Bi jn 为 mn 台 设 备 对 网 络 动态 的 修正 ，T(9)、。 系统 在 不 同 柄 态 下 的 稳定 性 ”总结 为 引 理 3。 具体 


7()、73(G) 和 7(s) 为 具有 合适 维 的 子 矩 阵 ; 万 为 证 明 见 附录 C。 


UE 


2 维 单位 矩阵 ; 0n，0. 分 别 为 n、 r 维 零 和 矩阵。 eA YN Ss peg a 
假设 设备 未 接 入 时 电网 是 稳定 的 。 根 据 Schur 引 理 2: 关注 模 态 下 广义 特征 抢 阵 满足 : 
定 理 P] ， det(T%a(s)=0 的 零点 由 det(Ti(s)- L.=Y" MX, =Yr, (Ss)+h ®E(s) (17) 


T2(s)T41(s)T3(3)) =0 决定 , 故 det(T”ha(3))=0 决定 了 式 中 ; 地 = 1, ® Yr (+(2/P)B ®¢(s) 
图 1 所 示 系 统 的 稳定 性 ， 其 B J | 
- ee 其 中 ，Yigreg(s) 为 相似 的 单 台 等 效 设备 动态 模型 ，k 
2 45) 为 第 上 个 关注 的 特征 模 态 。 
=det (Yn (s)+¥,(s))=0 当 设备 动态 模型 存在 差异 时 ， 可 构造 出 等 效 同 
构 系 统 , 其 矩阵 的 特征 子 空间 与 实际 系统 矩阵 的 特 


征 子 空间 近似 ， 由 此 得 到 广义 特征 矩阵 五 的 一 阶 
近似 。 具 体 地 ， 分 别 利用 办 ( 关 )、 加 (于 ) 构造 等 效 
同 构 系统 矩阵 Me， 满 足 : 

他 = (en ks ) + nom (5) (18) 
式 中 : yy (5s)=(U°/P)B', ®¢(s) 

Yr 1 (5)-1, ® (G(s)+ Be -p(s)) 


Cx (s) = 用 [diag ((U? /P )Yow (s))| XL 
> Pn (Uv? 2 Ti (s), 
Bes ss)=Y |(U°/P)(diag (Bo)+A 


Pn = Xa 


其 中 ，7Tiom x(5) 为 等 效 同 构 系统 和 矩阵，Xi= xxQP， 
Fi= yOL， 衣 、 天 的 列 问 量 分 别 为 (UP)B’ea @ 卫 
关于 妇 的 右 、 左 特征 向 量 ，xwz、yw 分 别 为 x、yi 的 
第 i 个 元 素 ; Gr (s)、Bceg$ 分 别 为 等 效 加 权 设 备 模 
型 与 等 效 加 权 电 容 容 纳 ;， pi 为 第 i 台 设 备 的 参与 因 
子 。 
每 个 等 效 同 构 系 统 中 设备 模型 相同 〈 均 为 等 效 
加 权 设 备 模型 )， 网 络 电容 也 均 为 等 效 加 权 电 容 ， 
其 余 网 络 动态 与 实际 系统 相同 。 将 等 效 加 权 设 备 与 
接地 电容 看 作 整 体 时 ， 根 据 引 理 1 和 引 理 2 可 知 ， 
1 的 一 组 特征 子 空间 仍 为 办) 、 只 (及 ) ， 且 它们 
是 二 特征 子 空间 的 零 阶 近似 ， 对 应 的 广义 特征 矩阵 
一 阶 近 似 分 别 表示 为 : 
L, = 了 4 MX, 
= GG, (s)+ Be “(Ss)+A, -CC(s) 


(19) 


» = 也 ((U°/P)B',® X 
让 由 0) ((V*/P)B', ®¢(s)) 
=yt (VU /P)B'x, 6(s) 
再 令 了, 一 阶 近似 的 逆 为 特征 子 系统 的 闭环 传 


递 函 数 , 那么 该 特征 子 系统 可 近似 实际 系统 在 对 应 
模 态 下 的 稳定 性 ， 即 定理 1， 具 体 证 明 见 附录 。 
定理 1: 将 实际 系统 矩阵 产 分 别 看 作 等 效 同 构 
系统 矩阵 Mi 的 摄 动 ，Mi 的 特征 子 空间 是 让 特征 子 
空间 的 零 阶 近似 ， 六 的 广义 特征 抢 阵 一 阶 近似 可 由 
式 计 算 。 
进一步 , 由 于 在 关注 模 态 下 应 的 广义 特征 矩阵 
近似 反映 了 实际 系统 的 稳定 性 ， 而 且 由 式 可 知 ， 广 
的 特征 子 系统 仅 为 两 输入 两 输出 系统 ， 因 而 基于 
该 特征 子 系 统 可 较 简单 地 实现 复杂 系统 的 分 析 和 
量化 。 
注 1: 当 设 备 动态 模型 同 构 ， 等 效 系统 矩阵 与 
实际 系统 矩阵 的 特征 子 空 = 间 恰 好 完全 相同 ， 特 征 子 
系统 稳定 性 与 原 系统 相同 ; 反之， 等 效 系统 矩阵 与 


(Bc )) ® nls) | XL 


实际 系统 矩阵 间 存 在 差异 , 广义 特征 矩阵 的 一 阶 近 
似 与 真实 广义 特征 矩阵 存在 偏差 ， 此 时 特征 子 系统 
近似 反映 了 系统 稳定 性 。 
3.3 ”特征 子 系统 的 物理 意义 

如 前 所 述 , 利 用 特征 子 空间 或 其 摄 动 可 分 别 求 
解 实际 系统 矩阵 的 关注 模 态 下 广义 特征 和 矩阵 或 划 
一 阶 近 似 ， 对 应 特征 子 系 统 可 近似 实际 系统 稳定 
性 ， 而 由 式 可 知 ， 由 于 等 效 网 络 模型 及 其 特征 子 
空间 的 结构 特征 ， 特 征 子 系统 物理 上 保留 了 设备 
与 网 络 的 动态 特性 ,恰好 为 低 维 单 馈 入 系统 ,因此 ， 
特征 子 系统 的 物理 意义 为 能 够 近似 实际 系统 稳定 
性 的 解 耦 子 系统 ， 如 图 2 所 示 。 

进一步 ， 由 式 与 图 2 可 知 ， 解 耦 子 系统 中 设 
备 动态 为 等 效 加 权 设备 动态 ,是 等 效 设备 动态 关于 
等 效 网 络 动态 的 加 权 平 均 ， 能 够 反映 等 效 设 备 动态 
的 整体 特性 ， 网络 动态 包含 5 等 效 线路 动态 及 设 
备 并 网 点 对 地 电容 动态 ， 其 中 ， 线 路 导 纳 为 不 同 模 
态 下 特征 子 空间 对 应 的 等 效 导 纳 ， 考 虑 了 设备 等 效 
导 纳 及 等 效 设备 实际 工 况 对 网 络 动态 的 修正 ， 能 够 
有 反映 系统 等 效 网 络 的 最 小 连接 强度 ; 对 地 电容 为 等 
效 加 权 电 容 ， 为 实际 系统 中 电容 关于 等 效 网 络 动态 
的 加 权 平 均 ， 反 映 实际 系统 电容 的 整体 特性 。 
恨 网 型 变 流 器 1 ”实际 复杂 系统 


加 


跟 网 型 ; 
了 / ntm 
-人 构 网 型 变 流 器 
必 “ -二 
一 ntmtr+1 外 部 电网 


特征 子 空间 | 和、 稳定 性 
近似 ”MM 等 价 


图 2 特征 子 系统 物理 意义 示意 图 


Fig.2 The schematic diagram of eigen-subsystem with its 
physical significance 


值得 注意 的 是 ,不同 模 态 对 应 的 特征 子 空间 及 
特征 子 系统 不 同 ， 对 应 磊 取 值 不 同 ， 而 在 实际 多 样 
化 设备 接 入 系统 中 ， 由 锁 相 环 带 来 的 小 扰动 同步 失 
稳 风 险 较 高 ， 尤 其 是 在 弱电 网 中 。 因 此 ， 经 常 关 注 
小 扰动 同步 稳定 下 的 特征 子 系统 (不 失 一 般 性 假设 
对 应 请 1)， 前 期 广义 短路 比 研究 主要 关注 该 问题 ， 
就 是 基于 该 特征 子 系统 近似 实际 系统 的 小 扰动 稳 
定性 。 其余 模 态 的 稳定 机 理 不 同 , 例如 电压 稳定 等 ， 
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对 应 大 可 能 不 同 ， 即 过 1， 但 仍 可 用 类 似 的 方法 近 
似 得 到 特征 子 系统 。 还 值得 一 提 的 是 ， 考 虑 更 复杂 
场景 时 ,如 设备 间 存 在 相 角 差 或 功率 因数 不 为 1 等 ， 
对 应 实际 系统 矩阵 、 特 征 子 空间 及 广义 特征 矩阵 为 
5 的 函数 ， 但 也 存在 对 应 的 特征 子 系统 。 针 对 其 余 
模 态 下 大 的 具体 取 值 、 复 杂 场 景 下 特征 子 系统 及 其 
应 用 ， 将 是 未 来 值得 深入 的 研究 工作 。 


4 ”基于 特征 子 系统 的 广义 短路 比 计算 


本 节 提 出 了 基于 特征 子 系统 的 广义 短路 比 通 
用 计算 原理 ， 并 提出 了 统一 计算 流程 。 此 外 ， 分 别 
给 出 了 几 种 典型 场景 下 的 指标 计算 方法 。 
4.1 广义 短路 比 的 通用 计算 原理 

如 前 所 述 ， 特 征 子 系统 的 物理 意义 是 实际 系 
统 的 解 而 子 系统 ， 它 由 等 效 加 权 设 备 动态 与 等 效 
线路 动态 构成 ， 可 近似 实际 系统 的 稳定 性 。 因 此 ， 
实际 系统 强度 取决 于 解 耦 子 系统 动态 ， 当 等 效 加 权 
设备 动态 一 定时 ， 实 际 系统 强度 取决 于 电网 广义 短 
路 比 ， 其 定义 为 第 大 个 关注 模 态 下 特征 子 空 间 对 应 
的 等 效 导 纳 ， 如 式 - 所 示 ; 进一步 ， 由 于 xt、yi 为 


(UWP)B’eq 关于 机 的 右 、 左 特征 同 量 ， 故 广义 短路 
比 等 于 为 特征 值 Mi， 其 表达 式 为 : 


2 
gSCRyr (U°/P)B'x = 马 汉 


Cg) 
ov a(oe 
其 中 ， 妈 (表示 求 第 大 个 特征 值 ， ca'x (.) 表 示 求 第 


个 奇异 值 ，Z Ye Be 。 

结合 上 文 分 析 与 式 可 知 ,， 不同 模 态 下 的 特征 子 
a 
而 前 期 广义 短路 比 主 要 关注 小 扰动 同步 稳定 模 态 
(=1)， 其 定义 就 是 特征 值 或 奇异 值 (和 1 或 oo )， 具 
体 见 文献 [21] 的 广义 短路 比 定义 。 其 余 重 要 的 特征 
子 系统 ， 即 k 关 1， 也 能 类 似 广 义 短路 比 定义 提出 指 
标 。 还 值得 特别 注意 的 是 ， 广 义 短路 比 既 可 以 通过 
和 矩阵 特征 值 定义 ， 也 能 通过 闭环 系统 的 特征 癌 量 得 
到 。 当 闭环 系统 的 特征 子 空间 与 开 环 系统 〈 或 电网 
等 效 导 纳 和 矩阵 ) 的 特征 子 空间 相同 时 ， 例 如 ， 同 构 
系统 ， 两 者 结论 一 致 ， 反之 ， 基 于 特征 值 的 定义 仍 
一 种 误差 较 小 且 便 于 实际 操作 的 简化 。 


广义 短路 比 越 小 ， 系 统 强度 越 弱 。 根 据 这 一 单调 特 
性 ， 能 够 确定 等 效 加 权 设 备 在 达到 临界 稳定 状态 或 
给 定性 能 要 求 所 需 的 最 小 短路 比 (简称 “等 效 加 权 
设备 临界 短路 比 SCRo”)。SCRo 可 反映 等 效 加 权 设 


了 结合 式 和 可 知 ， 当 等 效 加 权 设 备 动态 一 定时 ， 
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备 动态 对 电网 强度 的 耐 受 能 力 ， 恰恰 就 是 广义 短路 
比 临界 值 CgSCRo， 表 示 为 : 


SCR, = arg {det(G, (ss )+ gSCR-G(s,))=0) (21) 
其 中 ，arg{} 表 示 求 方程 的 根 ; sz 为 等 效 加 权 设 备 在 


给 守 中 


生 能 要 求 下 解 得 系统 主导 特征 根 。 

结合 上 述 过 程 ， 基 于 特征 子 系统 的 广义 短路 比 
计算 原理 如 图 3 所 示 。 有 具体 地 ， 利 用 计算 简单 且 元 
素 均 为 标量 的 等 效 时 纳 怎 阵 特征 子 空 间 , 能 够 求解 
实际 系统 矩阵 的 广义 特征 矩阵 一 阶 近似 ,对 应 特征 
子 系统 能 够 近似 实际 系统 稳定 性 ; 特征 子 系统 保留 
了 设备 与 网 络 的 动态 特性 并 实现 了 网 络 动态 与 设 
备 动态 的 分 离 ， 形 成 了 物理 意义 明确 、 融 合 电 网 广 
义 短 路 比 与 设备 临界 短路 比 的 系统 强度 量化 方法 。 
其 中 ，gSCR 仅 与 等 效 网 络 动态 有 关 ; 等 效 加 权 设 
备 临界 短路 比 SCRo 与 等 效 加 权 设 备 动态 有 关 。 


广义 短路 比 计算 原理 


实际 系统 矩阵 Tiei(s) 


基于 设备 物理 特 
性 、Schur 补 简化 
特征 子 空间 近似 


实际 系统 矩阵 下 = 7" (s) 


义 特征 矩阵 等 效 同 构 系统 矩阵 
介 近 似 Mo=Thom(s) 


特征 子 系统 : 解 耦 子 系统 
广义 特征 矩阵 
设备 电网 分 离 
等 效 加 权 设 备 
耐 受 能 力量 化 
SCRo 


有 网 强度 量化 
gSCR 


图 3 广义 短路 比 计算 原理 流程 

Fig.3 The process of gSCR calculation principle 
4.2 ”典型 场景 下 广义 短路 比 计算 方法 

基于 上 述 计 算 原 理 ， 本 节 给 出 三 种 典型 场景 中 
实际 系统 矩阵 的 具体 形式 ， 由 此 实现 广义 短路 比 的 
计算 。 

场景 1: 设备 在 非 额 定 运行 点 时 的 系统 广义 短 
路 比 。 考 虑 图 1 所 示 系 统 仅 包 含 n 台 采 me 
策略 、 参 数 的 新 能 源 变 流 器 ，m=0， 其 中 ， 新 能 
变 流 器 在 不 同 的 非 额定 运行 点 ， 和 输出 功 功率 及 端 
口 电 压 分 别 为 0.2~1p.u. 和 0.9~1.1p.u. 不 等 2g 。 
场景 2: 设备 功率 反 向 的 广义 短路 比 。 根 据 设 
备 的 物理 特性 与 系统 需求 ， 部 分 设备 可 反 向 从 电网 
吸收 功率 ， 如 储 能 系统 与 柔性 直流 输电 等 。 不 失 
般 性 地 ， 以 常规 储 能 系统 为 例 ， 考 虑 图 1 所 示 系 统 
包含 n 台 采 用 不 同 控制 策略 ,参数 的 新 能 源 变 流 嚣 ， 
m 台 常 规 储 能 系统 。 其 中 , 设备 均 运 行 在 额定 工 况 ; 


新 能 源 变 流 器 主要 向 电网 发 出 有 功 功率 ， 储 能 变 流 为 直观 展示 上 述 3 种 典型 场景 下 实际 系统 矩阵 
器 主要 向 电网 发 出 (放电 ) 或 吸收 有 功 功 率 (充电 )。 的 构成 , 表 1 分 别 给 出 不 同 场景 下 设备 动态 Y”1sr(s)、 

场景 3: 考虑 构 网 型 设备 的 广义 短路 比 。 考 虑 网 络 等 效 导 纳 和 矩阵 和 LP2/ 忆 等。 进一步， 结合 上 述 
1 所 示 系 统 包含 2 台 采 用 不 同 控制 策略 、 参 数 的 计算 原理 ， 即 可 分 别 求解 3 种 典型 场景 下 实际 系统 
新 能 源 变 流 器 ，m 台 构 网 型 储 能 系统 或 调 相机 。 其 和 矩阵 的 广义 特征 秆 阵 近 似 ， 得 到 特征 子 系统 ， 以 实 
中 ,设备 均 运行 在 额定 工 况 ; 调 相 机 和 外 特性 近似 现 广义 短路 比 及 设备 临界 短路 比 的 计算 。 此 外 ， 表 
为 电压 源 且 惯量 系数 较 大 的 构 网 型 变 流 器 能 够 为 1 还 给 出 了 上 述 场景 对 系统 强度 的 作用 机 理 ， 详 细 
0 在 锁 相 环 主导 的 小 扰动 同 。” 讨论 可 分 别 参 考 文献 [9][11][12][21][22]。 
步 稳定 问题 中 能 够 等 效 为 含 内 电抗 的 电压 源 。 

表 1 典型 场景 下 实际 系统 矩阵 构成 及 其 对 系统 强度 的 作用 机 理 
Tab.1 The elements required for calculation of system strength and their physical significance under typical 
scenarios 


实际 系统 矩阵 构成 ( 见 式 ) 
yu: BitBi BR 降 wan(9) 端 口 的 网 | 
人 FomaG) 0 We 
B. 设备 等 效 导 纳 Bi jpg 
非 额定 运行 点 diag( Yrsci(s)) 包含 所 有 设备 0 diag(U2/P) | 修正 等 效 网 络 动态 , 并 由 此 间接 修 
功率 反 向 diag(Yrsci(s) Yrcs(s)) | 及 中 间 节 点 的 Onurm+ SB 正 等 效 加 权 设 备 动态 , 从 而 分 别 修 
构 网 型 设备 diag( Yrsci(s)) 网 络 导 纳 矩阵 | diag(0, diag(B ap) 0;) Sp! 下 gSCR、SCRo 
表格 中 符号 及 公式 说 明 : 场景 1 中 0 系统 后 通过 测试 方法 获得 。 
i 台新 能 源 变 流 器 导 纳 传递 函数 和 矩阵， 二 1,.….,n。 场 5) 根据 式 -分 析 系 统 强度 及 稳定 裕 度 
景 2 中 Yrcy(s) 为 第 j 台 储 和 g 变 流 句 的 导 纳 传闻 男 数 值得 注意 的 是 ， 当 系统 为 黑箱 模型 时 ， 等 效 网 
矩阵 ;Saim =diag (Se So)，Sw 为 第 ;台新 络 动态 与 实际 工 况 矩阵 可 通过 扰动 辩 识 方法 gj 与 
能 源 变 流 器 的 额定 容量 , Sp0+j 为 第 j 台 储 能 变 流 器 。 ”PMU 等 量 测 数据 获得 
向 电网 发 出 或 从 电网 吸收 有 功 功率 的 额定 值 ， 以 发 综 上 所 述 , 通 过 引入 特征 子 系统 的 广义 短路 比 
出 有 功 功率 为 正 。 场 景 3 中 ，BL an 为 第 j 台 构 网 。 计算 方法 ， 统 一 了 多 种 场景 下 广义 短路 比 及 临界 
型 变 流 器 或 调 相 机 的 等 效 导 纳 ， 产 1 7， 下 的 计算 ， 具 有 优势 为 : 1) 保留 物理 特性 。 特 征 


BL Badiag(SvswjBp 或 1/xg”)，xg” 为 第 j 合 调 相机 子 系统 保留 了 设备 与 网 络 的 动态 特性 ; 2) 方法 通 
的 等 效 内 电抗 ， 即 次 暂 态 电抗 ，1/SvswBo 为 第 j 台 用 且 拓 展 性 好 。 不 仅 适 用 于 异 构 多 馈 入 系统 的 复 
构 网 型 变 流 器 的 等 效 内 电抗 ，Svsw 为 第 j 台 构 网 型 ” 杂 场 景 ， 而 且 有 潜力 推广 应 用 到 除 小 扰动 同步 稳 
变 流 器 的 额定 容量 ，B 为 单位 容量 下 构 网 型 变 流 器 。 ” 定 和 静态 电压 稳定 外 的 其 它 稳定 问题 。 
的 等 效 时 纳 。 区 区 

CE 5 算 例 分 析 及 验证 
人 首先 基于 4 馈 入 系统 验证 特征 子 系统 求解 方法 


es 


六 和 玫 比 的 。 及 其 在 描述 原 系统 稳定 性 的 有 效 性 ， 其 次 ， 措 建 革 
系统 强度 量化 实用 化 步骤 ， 具 体 如 下 : 实际 新 能 源 基地 的 电磁 暂 态 模型 ， 验 证 广义 短路 比 

1) 根据 设备 物理 特性 将 设备 动态 分 别 看 作 等 。 指标 在 新 型 电力 系统 强度 量化 时 的 有 效 性 ， 并 进 -一 
效 设备 动态 或 对 网 络 动态 的 修正 ; 步 分 析 设 备 在 非 额定 运行 点 、 功 率 反 向 和 包含 构 网 


2) 根据 网 络 参数 与 设备 特点 ， 基 于 式 得 到 等 ， 和 4 
型 设备 对 指标 计算 与 系统 强度 的 影响 规律 。 
效 网 络 动态 BL+Bz sg; 根据 实际 工 况 、 设 备 数 和 L . Mm i 本 
Bit Bi sx， 基于 式 获得 等 效 时 纳 矩 阵 (02/P)B oa 5.1 ”特征 子 系统 及 其 求解 方法 的 有 效 性 


®L; 国 El 

3) 根据 式 并 结合 所 需 分 析 模 态 ， 求 解 等 效 导 常规 储 色 
纳 和 矩阵 的 对 应 特征 值 或 奇异 值 ， 即 广义 短路 比 ; ME2 

4) 利用 等 效 导 纳 矩阵 关于 广义 短路 比 的 特征 风电 机 组 
子 空间 ， 根 据 式 求解 等 效 加 权 设 备 动态 ; 根据 式 并 5E3 
通过 解析 法 或 数值 计算 法 获得 SCRo。 有 具体 地 ， 当 本 
设备 / 场 站 的 控制 结构 及 其 参数 已 知 时 ，SCRo 可 解 构 疝 型 人 
析 计 算 ;， 若 具有 封装 模型 ， 可 构造 等 效 单 馈 入 仿真 图 44 馈 入 系统 


U 
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Fig.4 The topology of a four-converter System 

为 验证 特征 子 系统 求解 方法 及 其 在 描述 原 系 
统 稳定 性 的 有 效 性 ， 并 说 明 其 解 耦 子 系统 的 物理 意 
义 ,搭建 4 馈 入 系统 电磁 和 暂 态 仿真 模型 如 图 4 所 示 ， 
其 中 节点 1、4 为 常规 、 构 网 型 储 能 系统 接 入 节点 ， 
节点 2、3 为 新 能 源 变 流 器 接 入 节点 ， 设 备 控制 策 
略 、 参 数 及 网 络 参数 如 附 表 E1-E3 所 示 。 
具体 地 ， 考 虑 等 比例 增 大 附 表 E3 所 示 所 有 线 
路 电感 ， 比 例 参数 三 1、1.5、2、3、4， 其 中 ， 常 
规 储 能 变 流 器 从 电网 吸收 功率 ， 构 网 型 变 流 器 向 电 
网 发 出 功率 , 设备 实际 工 况 如 表 2 所 示 。 以 此 为 基 
础 ， 首 先 ， 分别 基于 特征 子 空间 摄 动 理论 求解 特征 
子 系统 ， 其 次 ， 对 比 原始 4 馈 入 系统 和 特征 子 系统 
的 主导 特征 根 ， 验 证 特征 子 系统 的 有 效 性 。 
首先 ， 以 算 例 1 为 例 ， 给 出 特征 子 系统 的 具体 
求解 过 程 。 构 建 4 馈 入 系统 的 闭环 传递 函数 逆 矩 阵 
并 基于 物理 特性 简化 ， 有 具体 步骤 为 考虑 构 网 型 变 流 
器 的 等 效 导 纳 对 网 络 的 修正 ， 进 而 基于 Schur 定理 
将 4 馈 入 系统 矩阵 降 阶 为 仅 包 含 等 效 设 备 母 线 的 矩 
阵 ， 并 进一步 用 (U2/P)8@ 荆 左 乘 算 阵 , 最 终 4 馈 入 系 
统 矩 阵 变换 为 : 

det (Ti (s))= det (YH (s)+7"(s))=0 (22) 

式 中 : yw(s)= (U0°/P)B, ®¢(s) 


, | U2 U2 U2 
了 mx (diag Gr (s) Ysci (s) Vses (9 
1 


已 五 

-55 14 14 
(2/P)Bw =| -2 7 2 

-4 -4 16 


式 中 ，Yysci (s)、Yyscz (s)、Ypcs(s) 分 别 为 新 能 源 变 
流 器 与 常规 储 能 变 流 器 的 导 纳 传递 函数 和 矩阵 ;U2/P 
=diag(UI2PI LU22%P LU32P3)。 

等 效 导 纳 和 矩阵 (LP2/ 疡 Be @ 五 有 三 个 二 重 特征 
值 4=5.19、 ?2=17.10、43=-53.62， 对 应 有 特征 子 空 
间 R( 兰 )、R(Y)G=1,2,3), 满足 下 式 , 其 中 六 =w@ 
7， 天 六 ® 7，Xi、 二 分别 为 (U2/P) Beq 关于 的 右 、 
左 特征 向 量 。 

Y”((U*/P)B', ®L,)X,=4 8 (23) 
则 基于 特征 子 空间 摄 动 理论 ， 利 用 R(X)、 
R( 瑟 可 将 4 馈 入 系统 矩阵 近似 解 耦 为 三 个 广义 特 
征 和 矩阵 ， 表 示 为 : 

LYT (SX =G(s)t4.6(s) 024) 

式 中 3 py 二 Xi Vy 


马 


2 2 


U U U 
Ci = pi PY (s)+ Pi2 Pp Ye (s)+ Pis pe (s) 


3 


其 中 ， 庆 1,2,3，xy、 屹 分别 为 x、yi 的 第 j 个 元 素 。 
表 2 算 例 的 比例 参数 及 实际 工 况 
Tab.2 The proportional parameters and actual 
operating conditions of examples 


算 例 | U P SSCR 

1 1 | 1.05、0.99、1.03 5.19 

2 1.5 1.04、0.97、1 3.39 

3 2 | 1.03、0.95、0.99 | -0.2、1、0.6 2.47 

4 1.01、0.95、0.98 1.67 

5 4 | 0.99、0.93、0.97 1.21 
四 馈 入 系统 矩阵 的 特征 函数 与 这 三 个 广义 特 


征 和 矩阵 的 特征 函数 近似 相同 ，》 为 等 效 导 纳 矩 阵 的 
最 小 正 特征 值 ， 则 对 应 天 为 实际 系统 的 解 耦 子 系统 


且 能 够 近似 实际 系统 稳定 性 ， 令 万 为 特征 子 系统 
的 闭环 传递 函数 , 理论 上 该 特征 子 系统 可 近似 实际 
系统 稳定 性 。 进 一 步 ， 图 5 给 出 算 例 中 原始 4 馈 
入 系统 和 特征 子 系统 的 主导 特征 根 ， 可 以 看 出 ， 随 
着 电网 强度 的 变化 ， 原 始 4 馈 入 系统 和 特征 子 系统 
的 特征 根基 本 保持 一 致 ， 因 此 验证 了 特征 子 系统 可 
近似 实际 系统 稳定 性 。 
O 〇 4 馈 入 系统 


X 特征 子 系统 


y 


gSCR=1.21 


-S$ -4 -3 -2 -1 0 1 2 
图 5 系统 等 效 前 后 主导 特征 根 对 比 图 
Fig.S Comparison of dominant roots before and after 


equivalence 


5.2 ”基于 广义 短路 比 的 系统 强度 量化 
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图 6 某 实际 新 能 源 基地 拓扑 图 
Fig.6 The topology of a practical renewable energy plant 


为 验证 广义 短路 比 指标 在 新 型 


电力 系统 强度 


量化 时 的 有 效 性 ， 搭 建 某 实 际 新 能 源 基 地 的 电磁 暂 
态 仿 真 模型 ， 如 图 6 所 示 。 该 系统 共 包 括 4 个 区 域 
k 27 个 风力 发 电机 组 、2 个 常规 储 能 系统 和 2 个 构 
网 型 储 能 系统 ， 其 中 风力 发 电机 组 或 储 能 系统 的 下 
标 数字 代表 其 内 部 采用 变 流 器 编号 ; 图 中 节点 
1~54 表示 设备 接 入 节点 ， 其 中 节点 43、54 为 构 网 
型 变 流 器 接 入 节点 ， 节 点 14、15、28、29 为 常规 
储 能 变 流 器 接 入 节点 ， 其 余 节 点 为 新 能 源 变 流 器 接 
入 节点 ;节点 55~92 表示 中 间 节 点 ， 节 点 93 为 
场 站 外 电网 等 值 节点 。 有 具体 地 ， 考 虑 如 下 四 个 算 例 
情况 ， 算 例 的 设备 控制 策略 、 参 数 及 网 络 参 数 均一 
致 ， 如 附 表 E1-E3 所 示 ; 设备 实际 工 况 如 附 表 E4 
所 示 。 

算 例 1: 系统 中 设备 均 运 行 在 实际 工 况 ， 即 输 
出 有 功 功 率 及 端口 电压 往往 偏离 额定 值 ， 其 中 ， 常 


规 储 能 从 电网 吸收 功率 ， 其 余 设备 向 电网 发 出 功率 ; 


算 例 2: 在 算 例 1 的 基础 上 将 设备 工 况 均 调整 
为 额定 工 况 ， 其 中 ， 和 常规 储 能 仍 从 电网 吸收 功率 ， 
其 余 设 备 向 电网 发 出 功率 ; 

算 例 3: 在 算 例 1 的 基础 上 将 常规 储 能 变 流 器 
的 工作 模式 改变 为 发 出 功率 ; 

算 例 4: 在 算 例 1 的 基础 上 切除 构 网 型 变 流 器 ， 
即 系统 中 构 网 型 变 流 器 接 入 节点 悬空 。 


图 7 设备 端口 电压 的 时 域 波 

Fig.7 Time domain waveform of terminal voltages of 

devices in a practical renewable energy plant under 
different operating conditions 


首先 ， 按 照 上 述 指标 计算 流程 分 别 解析 计算 四 
个 算 例 中 系统 的 广义 短路 比 ; 其 次 ， 基 于 式 并 通过 
搭建 等 效 加 权 单 机 仿真 或 解析 方法 得 到 四 个 算 例 
中 系统 的 SCRo; 最 后 ， 基 于 式 ， 通 过 比较 gSCR 与 
SCRo 判断 系统 稳定 性 ， 并 进一步 基于 式 计 算 p%， 
以 评估 系统 强度 。 表 3 给 出 指标 值 及 评估 结果 。 可 
以 看 出 ， 算 例 1、3、4 中 gSCR 均 大 于 SCRo， 系 统 
满足 小 扰动 稳定 的 最 低 要 求 ， 其 中 ， 算 例 1 中 
po%>20% ， 满 足 小 扰动 稳定 裕 度 需求 ; 算 例 2 中 
gSCR=2.39<2.63=SCRo， 系 统 失 稳 。 

相应 的 , 在 仿真 时 间 T=0.3s 时 , 在 无 穷 大 母线 
处 (节点 93) 施 加 了 5% 端 电压 又 升 并 快速 清除 ， 售 
例 1、2 中 设备 端口 电压 波形 分 别 如 图 7(a)(b) 所 示 ， 
其 中 , 算 例 1 中 的 设备 端口 电压 振荡 快速 收敛 稳定 ; 
算 例 2 中 的 设备 端口 电压 振荡 发 散 ， 系 统 失 稳 。 上 
述 时 域 仿真 趋势 与 基于 广义 短路 比 的 系统 强度 评 
佑 结果 趋势 一 致 ， 故 广义 短路 比 在 系统 强度 量化 评 


NY 


表 3 算 例 的 指标 数值 
Tab.3 The values of indexes under four scenarios 
算 例 gSCR SCRo [00 
1 3.28 2.67 22.8% 
2 2.39 2.63 -9.1% 
3 3.03 2.68 13.1% 
4 3.02 2.69 12.3% 
1.15 
| 
S 
:1.05 
局 
1 


估 中 的 有 效 性 得 到 验证 。 
此 外 ， 图 7(c) 给 出 算 例 1 中 系统 在 gSCR= 
SCRo 时 设备 端口 电压 的 等 幅 振荡 波形 ， 验 证 了 
gSCR= SCRo 在 描述 系统 稳定 边界 的 有 效 性 ， 此 时 
系统 网 络 参数 如 附 表 E3 所 示 。 
5.3 ”广义 短路 比 的 影响 因素 
5.3.1 ” 非 额定 运行 点 对 指标 计算 影响 
为 揭示 非 额定 运行 点 对 指标 计算 及 系统 强度 
9 影响 规律 ,考虑 对 比分 析 设 备 在 实际 工 况 ( 算 例 1) 
与 额定 工 况 的 场景 ( 算 例 2)。 
首先 ， 在 指标 计算 过 程 中 ， 算 例 1、2 的 区 别 
仅 在 于 和 矩阵 U2P 的 不 同 ， 其 中 ， 额 定 工 况 下 LPZ/P=- 
SBs21，SBs2 为 52 X52 的 对 角 和 矩阵 ， 元 素 为 设备 的 
额定 容量 , 以 发 出 有 功 功 率 为 正 ; 实际 工 况 下 UP= 
diag(U2/Pi)s2，diag(U2/Pi)sz 为 52X52 的 对 角 和 矩阵 。 
其 次 ， 比 较 算 例 1、2 中 广义 短路 比 及 SCRo 的 
变化 。 可 以 看 出 ， 两 算 例 中 广义 短路 比 与 SCRo 均 
不 同 ， 验 证 了 非 额 定 运 行 点 会 同时 改变 广义 短路 比 
及 SCRo， 从 而 影响 系统 强度 ， 说 明了 考虑 实际 工 


i 


第 于 刘 居 吐 等 ， 基 于 特征 子 系统 的 广义 短路 比 计算 原理 11 
况 对 系统 强度 影响 的 必要 性 。 一 步 拓展 到 高 频 模 态 ， 故 它 具 有 广阔 的 应 用 前 景 
5.3.2 ”功率 反 向 对 指标 计算 影响 参考 文献 
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理 


流 右 接 入 ( 算 例 


的 gSCR 得 到 一 定 提升 ， 而 SCRo 改变 较 小 ， 验 
了 储 能 变 流 器 工作 于 充电 模式 会 提升 广义 短路 
， 进 一 步 ， 通 过 改变 等 效 网 络 动态 中 算 阵 UW/P 
间接 改变 等 效 加 权 设 备 动 态 ， 从 而 修正 SCR。。 
影响 
为 分 析 构 网 型 设备 的 影响 规律 ， 考 虑 构 网 型 变 
1) 与 不 接 入 的 场景 ( 算 例 和 。 

首先 ， 根 据 构 网 型 设备 的 物理 特性 可 知 ， 在 指 
计算 过 程 中 ， 算 例 1、4 的 主要 区 别 在 于 设备 对 
等 效 网 络 动态 的 修正 Bi 75R 不 同 ， 其 中 ， 考 虑 构 网 
型 变 流 器 接 入 时 ，Bz zsr 中 与 构 网 型 变 流 器 接 入 节 
对 应 的 对 角 元 为 构 网 型 变 流 器 的 等 效 导 纳 ， 会 修 
FE 网络 动态 ， 而 不 接 入 时 Bz zsr 中 元 素 均 为 0。 
其 次 ， 比 较 算 例 1、4 中 广义 短路 比 及 SCRo0。 
以 看 出 ， 接 入 构 网 型 变 流 器 后 系统 广义 短路 比 得 
一 定 提升 ， 而 SCRo 改变 较 小 ， 说 明 构 网 型 变 流 
的 接 入 提升 了 交流 电网 强度 。 


结论 与 展望 


1) 给 出 了 多 馈 入 系统 的 特征 子 系统 定义 、 物 
意义 及 计算 方法 ， 在 此 基础 上 提出 了 广义 短路 比 


及 
络 


太 


其 临界 值 的 通用 计算 方法 。 特 征 子 系统 保留 了 网 
与 设备 的 动态 特性 并 实现 了 网 络 动态 与 设备 动 


态 的 分 离 ， 基 于 广义 短路 比 的 系统 强度 量化 方法 物 


理 意义 清晰 、 实 际 应 用 时 方便 。 

2) 基于 特征 子 空间 摄 动 理论 ， 分 析 了 广义 短 
路 比 计算 的 误差 来 源 , 并 针对 跟 网 型 设备 在 非 额定 
运行 点 、 有 功 功率 反 向 和 系统 包含 构 网 型 设备 等 
型 实际 应 用 场景 ， 给 出 了 广义 短路 比 及 其 临界 
直 的 计算 方法 。 

新 型 电力 系统 有 多 个 特征 子 系统 ， 其 物理 意义 
和 对 应 的 稳定 问题 不 相同 。 虽 然 本 文 聚焦 低频 模 态 


的 


特征 子 系统 ,但 基于 特征 子 系统 区 


思路 也 可 以 进 
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附录 A 引 理 1 的 举例 
以 两 维和 矩阵 为 例 说 明 引 理 1，M1o 与 万 分 别 为 : 


1 2 | 2 . 
M,= 2 ® 1,, M= pe ® ,+diag(s,) (25) 


其 中 ，s;E R， 数 值 较 小 ， 二 1,2,3,4。 
根据 式 ， 利 用 Wo 的 两 个 特征 子 空间 可 求解 六 
广义 特征 矩阵 的 一 阶 近似 : 


下 


L 人 Y” MX, 
0.75e +0.24e, (26) 
=—0.37@®L+ 
0.75e,+0.24e, 
L 人 YMX, 
0.24e +0.766, (27) 
=5.37®L+ 
0.24e, +0.76é, 


式 中 ; X=[-0.82 057] 8@L,%=[-0.42 -091] @P 
了 =[-0.92 0.42]@ L,Y =[-0.57 -0.84]® 

中 ,入 |、 各 中 的 列 向 量 为 Mo 关于 二 重 特征 值 -0.37 

与 5.37 的 右 特征 向 量 ,， 五 、 钱 中 的 列 向 量 为 对 应 左 

特征 向 量 ， 可 分 别 张 成 特征 子 空间 。 

如 果 e1、s2、s3、&4 分 别 为 0.4、0.3、0.2、0.1， 

义 特征 矩阵 及 其 一 阶 近 似 的 结果 分 别 为 : 


| 由 


y 


忌 ， 


第 刘 晨 咯 等 ; 基于 特征 
加 pal os) 
有 bs 区 BS | (29) 
—0.121 $5.518 
附录 B 式 的 详细 推导 
高 维 实际 系统 矩阵 为 : 
| 了 s) Tls 
| G0 
式 中 Z(s)=B@6(s),T(s)=B OC(s) 
T(s)= Yaa(s )+B el(s )+diag(B. ; ) Bls) 
T,(s)= B,®¢el(s )+ diag (Be ,)® 1(s) 
其 中 ,i=1,...yn, 太 ntl,..., ntmtro 
则 Ti(s)- T2(s)T41(s)T3(s) 可 表示 为 : 
7 (5)-1(s) 77 (s)L,() 
= Ys(s)+ B®E(s)+diag (Bs) Ou(s) 31) 


-(B,B; ®1,)K(s)(B, 6(s)) 


式 中 : K(s)= (1 + diag (Bo) Br ® p(s)6 (5)) 
根据 Neumann 级 数 性 质 可 知 ， 对 于 矩阵 4(s)， 
若 其 谱 半 径 满 足 pCd(J)<1，CH A(3))1 守 7 -4(9)， 
考虑 线路 接地 电容 较 小 ， 则 K(s) 可 转化 为 : 
K(s)~ Tn) -diag (Bo) Bi Bu(s)e(s) (32) 
将 式 代 入 式 可 得 : 
T'(s)-T(s) TH (s)T(s) G 
= Yar(s)+ B's @6(s)+(diag (Bc)+A(Be))® p(s) 


Bo )Bi'B, 


| 
ml 


式 中 : A(B8.)= 
附录 C 引 理 2 的 证 明 

根据 特征 值 定义 可 知 (UWP)B’eq 8 有 nn 个 二 
重 特征 值 4(i=1,.….,n); 对 应 (UWP)B’”en BD 有 一 组 特 
征 子 空间 RR( 关 )、R(Y)， 其 中 , %=xi@ ，Y8= 
J 站 8 卫 ，XX、 了 i; 列 癌 量 分 别 为 等 效 导 纳 和 矩阵 关于 扩 
的 右 、 左 特征 向 量 。 则 利用 吕 兰 ) 和 中 办 可 求解 
(DWP)B’eaq 8 也 的 n 个 广义 特征 矩阵 : 

L =Y"((U°/P)B', BL )X,=48I (34) 

当 实 际 系统 矩阵 让 中 等 效 设备 模型 相似 ， 均 、 

玉生 满足 如 下 等 式 : 
MX =XL,Y M=LY 

式 中 : =ou()+4@56(s) 
再 结合 式 与 定义 12 可 知 ，R( 半 )、R(FP) 和 了 工 
分 别 为 实际 系统 矩阵 的 特征 子 空间 及 广义 特征 先 
阵 。 进 一 步 ， 令 (Xi ... 加 )、 天 (Fi ... 了 ,)， 利 用 
XX、 了 可 将 实际 系统 矩阵 转换 为 : 


B,B, 1diag ( 


(35) 


子 系统 的 广义 短路 比 计算 原理 


下 = 了 ”MX = diag (L,) 


13 
(36) 


由 于 惟 了 的 列 向 量 均线 性 无 关 ， 故 det( 万 )-0 


Sdet(£ LZ)= 0 等 价 ， 再 


为 : 


det(Z)= oTlaet(L 中)= 0 
] R(X)、R(Y) 可 将 实际 系统 条 


阵 分 别 解 砖 为 n 个 广义 特征 矩阵 了， 进一步 ， 结 合 
式 可 知 ， 广 义 特征 矩阵 保留 了 等 
所 以 式 可 理解 为 将 实际 系统 解 看 为 n 个 短 
馈 入 系统 。 因此 ， 利 


mm(Z ) 可 求解 五 ， 对 应 的 特征 子 系统 


上 式 表明 , 利 


态 特 性 ， 
路 比 不 同 的 等 效 单 


由 式 可 知 ， det(£)= 0 可 转化 


(37) 


效 单 馈 入 系统 的 动 


1 R(X,)、 


为 特 


征 子 系统 的 闭环 传递 函数 ) 可 近似 实际 系统 在 对 应 
模 态 下 的 稳定 性 。 证 毕 


附录 D 定理 1 的 证 明 
证 明 : 利 


1 RX). RK 
构 系 统 矩 阵 Mi 的 广义 特征 和 
LY MX,=G(s)+ 


中) 可 求解 分 别 等 效 同 


E 阵 ， 满足 : 
Bia HS) 6(s) (38) 
结合 式 、 可 知 ， 等 效 同 构 系 统 的 广义 特征 矩阵 


是 实际 系统 的 广义 特征 矩阵 一 阶 近 似 。 根据 引 理 1， 
将 实际 系统 矩阵 看 作 等 效 同 构 系统 矩阵 的 摄 动 ， 故 
利用 只 总 )、R(¥ ) 即 可 得 广义 特征 年 阵 的 一 阶 近 
似 。 证 毕 。 
附录 了 
表 El 新 能 源 及 储 能 变 流 器 控制 参数 
Table El Control parameters of grid-following converters 
内 容 新 能 源 变 流 器 储 能 变 流 器 
交流 系 统 额定 容 1500kVA, 620V 
量 、 电 压 
人 0.68MVA(14,15) 
设备 额定 容量 1.5MVA O70 NYA a 70) 
直流 侧 电容 、 电 压 0.1089p.u., 1.1kV 0.1089p.u., 1.5kV 
外 环 PI 0.5,40( 直 流 电压 ) 0.5,20( 功 率 ) 


滤波 电感 、 电 容 、 


0.05、0.0006、0.002 


0.2,60 


0.05、0.0006、0 


0.5,15(14,15) 
0.3,10(28,29) 


四 馈 入 系统 : 14,10400 
某 实 际 新 能 源 场 站 : 


锁 相 环 PI 


10,10500(1-13); 21,10300 


40,24000 


(16-27); 26,10400 (30-38); 
16,10200(39-42,44-48); 
31,10200(49-53) 


电压 前 馈 滤波 时 间 


常数 0.00001 
表 E2 构 网 型 变 流 器 控制 参数 


Table E3 Control parameters of grid-forming converter 


内 容 


数值 


设备 额定 容量 


Bond 
也 


电压 外 坏 PI 参数 
滤波 电感 、 滤 波 电 容 、 
电流 内 环 PI 参数 

构 网 型 变 流 器 惯量 系数 、 阻 尼 系 数 
电压 前 馈 滤波 时 间 常 数 


电阻 


1.S5MVA 
0.8,12 


0.05、0.0005、0.02 


6,10 
2,20 
0.001 


表 E3 仿真 电网 参数 


Table E3 Parameters of Simulations 


I 馈 入 系统 电网 参数 (p.u.) 


X1.s=0.20 X26=0.40 X37=0.60 X12=0.40 X23=0.40 
X13=0.40 X34=0.20 

某 实际 新 能 源 基 地 电网 参数 (X103p.u.) 
Xsss6=0.24 Xs6.86=1.00 Xs7.86=2.40 Xss,86=3.30 Xso.86=2.90 
X6086=2.00 X61.86=1.70 X62.86=0.03 X63.87=1.10 X64.87=3.10 
X65.87=2.30 X66.87=1.90 X67.87=1.40 Xeg,87=0.35 X69.87=0.16 
X70.87=0.02 X71.88=2.50 X72.88=2.20 X73,.88=1.50 X74.88=1.30 
X75,88=0.81 X76.88=1.30 X77.88=0.81 X78.90=0.19 X79.90=2.20 
X80.89=0.19 X8189=2.20 X8291=1.50 X8391=1.30 X8491=1.50 
Xs8s01=1.30 X86.87=1.00 X87.92=3.10 Xsgg,92=4.40 Xsg9.90=7.00 
Xo0.02=7.80 Xo1.0=0.87 Xo2.93=6.90 变压器 电抗 : 64.0 


算 例 1 在 gSCR=SCRo 时 的 电网 所 有 参数 是 上 述 电 网 参数 的 1.1 倍 


表 E4 某 实际 新 能 源 基地 4 个 算 例 中 的 设备 实际 工 况 


Table E4 The operating conditions of devices in four examples 
of Practical Renewable Energy Plant 
算 例 运行 工 况 
Pp 1-13:0.7;16-27:0.8;30-38:0.6;39-42:0.75; 44-53:0.9;14-15: 
1 -0.45;28-29: -0.48;43,54:1.00; 
7 1:13, 16-27, 30-41, 44-48 :1.01;14-15, 28-29:1.03;42,43,54:1.00; 
49-53:1.02 
2 1.00 
U 1.00 
3 2 14-15: 0.45;28-29: 0.48; 其 余 见 算 例 1 
U 14-15, 28-29:1.01， 其 余 见 算 例 1 
了 P 见 算 例 1 
U 见 算 例 1 
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